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ABSTRACT 

Using Herschel PACS and SPIRE observations of Lockman Hole-North and GOODS- 
N as part of the HcrMES project, we explore the far-IR properties of a sample of 
mid-IR selected starburst dominated ultra-luminous infrared galaxies (ULIRGs) at 
z ~ 2. The selection of the sample is based on the detection of the stellar bump 
that appears in the SED of star-forming galaxies at 1.6/xm. We derive robust esti- 
mates of infrared luminosities (£ir) and dust temperatures (Td) of the population 
and find that while the luminosities in our sample span less than an order of magni- 
tude (12.24 < log(LiB,/L Q ) ^ 12.94), they cover a wide range of dust temperatures 
(25 ^ Td ^ 62 K). Galaxies in our sample range from those that are as cold as high- 
z sub-millimeter galaxies (SMGs) to those that are as warm as optically faint radio 
galaxies (OFRGs) and local ULIRGs. Nevertheless, our sample has median Td=42.3 K, 
filling the gap between SMGs and OFRGs, bridging the two populations. We demon- 
strate that a significant fraction of our sample would be missed from ground based 
(sub)mm surveys (850-1200/im) showing that the latter introduce a bias towards the 
detection of colder sources. We conclude that Herschel observations, confirm the exis- 
tence of high-z ULIRGs warmer than SMGs, show that the mid-IR selection of high-z 
ULIRGs is not Td-dependent, reveal a large dispersion in T^ of high-z ULIRGs, and 
provide the means to characterize the bulk of the ULIRG population, free from selec- 
tion biases introduced by ground based (sub)mm surveys. 



1 INTRODUCTION 



One of the most successful methods for selecting high-z 
ultra-luminous infrared galaxies (ULIRGs: I/8-1000 ^m > 
10 12 Lq) is their direct far-IR detection via ground based 
(sub)millimeter surveys (e.g. Barger et al. 1998, Hughes et 
al. 1998, Mortier et al. 2005, Pope et al. 2006, Austermann 
et al. 2010). This technique has revealed the population of 
the so called submillimetre galaxies (SMGs), that represent 
a significant class of high-z ULIRGs. Attempts to charac- 
terize their dust temperature (Td) show that these galax- 
ies are colder when compared to local ULIRGs (e.g Chap- 
man et al. 2005), suggesting that in general high-z ULIRGs 
tend to have lower dust temperatures. However, the sub- 
millimetre technique introduces a bias towards the selection 
of ULIRGs with lower dust temperatures while it misses 
warmer ULIRGs. First observational evidence of a missing 
population of high-redshift dusty star-forming galaxies with 
hotter dust has been given by Chapman et al. (2004) using 
a selection of radio-detected but sub-mm-faint galaxies with 
UV spectra consistent with high-redshift starbursts. These 
optically faint-radio galaxies (OFRGs) share similar prop- 
erties with SMGs (e.g. stellar mass, SFR) but some have 
considerably higher dust temperatures (Casey et al. 2009, 
Magnelli et al. 2010). 

Another technique that has been proven to pick high-z 
starburst dominated ULIRGs efficiently, is based on mid-IR 
color selection. This technique relies on the detection of the 
rest-frame 1.6/im bump in the SED of star-forming galax- 
ies, produced by thermal emission from late-type stars and 
enhanced by an apparent emission feature due to H- ions 
in the atmospheres of giant stars (Simpson & Eisenhardt 
1999; Sawicki 2002). The advent of Spitzer allowed the de- 
tection of this feature in z~2 galaxies and subsequent IRS 
spectroscopy has demonstrated the efficiency of the method 
to select star-burst dominated ULIRGs in a redshift range 
of 1.5 < z < 2.5 (e.g. Farrah et al 2008, Weedman et al 

2008, Huang et al. 2009). Further studies (Lonsdale et al. 

2009, Fiolet et al. 2009, Kovacs et al. 2010), indicate that 
only 40% of their sample is made up of bright mm sources 



and thus belong to the class of SMGs, while most of the rest 
have lower Si. 2mm fluxes. 

The above suggests that a considerable fraction of mid- 
IR selected high-z ULIRGs are missed by ground based 
(sub)-mm surveys. While this could naturally be explained 
if their dust temperature (for a given luminosity) is higher 
than that of the SMGs, this is not yet clear, as up to now 
the study of their far-IR properties is restricted to objects 
with ground (sub)mm detection or to the most luminous ex- 
amples of the population with the highly confused BLAST 
beam (e.g Ivison et al. 2010, Dunlop et al. 2009). Hence, a 
far-IR study of the population, free of the selection bias in- 
troduced by the ground based sub-nim detection is required. 
Furthermore, detailed study of high-z ULIRGs is essential 
as up to now theoretical models fail to account for the in- 
ferred luminosities, star formation rates and number counts 
(Baugh et al. 2005, Dave et al. 2010). 

In this study, we use observations of Lockman-North 
(LHN) and GOODS-N fields obtained by the Herschel Space 
Observatory (Pilbratt et al. 2010) as part of the Herschel 
Multi-tiered Extragalactic Survey (HerMES, Olivier et al. 
2010), to investigate the FIR properties of a sample of mid- 
IR selected ULIRGs at z~2 (IRAC peakers). Taking ad- 
vantage of both the Photodetector Array Camera (PACS, 
Poglisch et al. 2010) and the Spectral and Photometric 
Imaging Receiver (SPIRE, Griffin et al. 2010) data that 
probe the peak of the SED of galaxies at this redshift, we 
derive robust dust temperature measurements for the bulk 
of the population and compare our sample to that of other 
high-z ULIRGs. Throughout this paper we assume fi m =0.3, 
H =71km sec -1 Mpc -1 and fi A =0.7 



2 SAMPLE SELECTION AND HERSCHEL 
OBSERVATIONS 

To select our sample we adopt the IRAC color criteria in- 
troduced by Huang e al. (2009). In particular we search for 
galaxies in LHN that satisfy the following IRAC color cri- 
teria: 0.05<[3.6]-[4.5]<0.4 and -0.7<[3.6]-[8.0]<0.5, have 
f24>0.2mjy and r veg a>23.0 to avoid low redshift interlopers. 
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Figure 1. Rest— frame SEDs and derivation of the far-IR properties for two ULIRGs in our sample. Solid orange line shows the best fit 
template up to observed 8/xm as derived by LEPHARE photo-z code. Solid black line shows the best fit CE01 model while the blue line 
depicts the best-fit modified black body (with /3=1.5), used to derive T<j estimates. The vertical dotted line indicates the wavelength cut, 
below which photometric data where not considered in the modified black body fit. Red circles denote it Herschel data. 



These criteria probe the 1.6/^m stellar bump while the red 
color cuts ensure the rejection of power-law AGNs. Subse- 
quent IRS spectroscopy of ULIRGs selected with the above 
method, has shown that this selection picks ULIRGs at very- 
narrow redshift range 1.7<z<2.3 with strong PAH features, 
indicative of intensive star-formation (Huang et al. 2009). 

We use the Spitzer Wide-Area Infrared Extra- 
galactic Survey (SWIRE) multi-wavelength catalog 
(U,G,R,I,z,J,H,K+IRAC+MIPS) (Surace et al. 2005) over 
the 0.25(sq deg) of LHN covered by PACS and SPIRE, 
and we identify 32 objects that meet our criteria. We then 
match the sample with the joined Herschel PACS (100- 
and 160^m and SPIRE (250- 350- and 500^m) XID catalog 
(Roseboom et al. 2010). For the PACS data, where source 
confusion is less severe, we consider fluxes based on blind 
source extraction. For the SPIRE data we adopt the fluxes 
derived with source extraction based on 24/im priors, we 
reject candidates with neighboring (d<20") 24/xm sources 
whose /24 is > 50% than that of our object. Finally, we 
require at least one PACS and SPIRE detection (3a). 

In the resulting LHN sample there are 25 candidate 
ULIRGs at z~2, five of which have spectroscopic redshifts 
(Fiolet et al. 2010, in prep). For the rest, we derive photo- 
metric redshifts, using the LEPHARE photo-z code (Ilbert 
et al. 2009). Namely, we fit the SED of the galaxies up to 
8.0/im with a wide range of template SEDs and consider dust 
attenuation that follows the prescription of Caltzetti et al. 
(2000) . For each object we adopt the redshift corresponding 
to the minimum \ 2 value of the fit. Two examples of the best 
fit template for two galaxies in our sample are depicted in 
Figure 1. The uncertainty of the photometric redshifts was 
derived based on the redshift probability distribution func- 
tion ( PDF(z) ) and we choose to exclude from our analysis 
candidates with multiple solutions or uncertainties larger 
than Az — 0.5. Furthermore, a comparison of the derived 
photometric redshifts with the spectroscopic redshifts that 
is available for 5 sources, yields a very good agreement be- 
tween the two values Az-=(z photo — z spcc )/(l+z spcc ) < 0.1. 
The final LHN sample consists of 18 ULIRGs with median 
z=1.98, and range 1.5<z<3.0. Finally, we match this sample 
to radio VLA 1.4GHz catalog of LHN (Owen et al. 2008). 



To increase the size of our sample, we perform the 
same procedure in the GOODS-N field. Using the multi- 
wavelength catalog and the SPIRE data (PACS data for 
GOODS-N are not available for this study), we indentify 
candidate ULIRGs with at least two detections (3a) at 
SPIRE bands. To exclude sources with strong AGN ac- 
tivity candidates with X-ray detection (I/x[0.5-8.0keV] > 
3xl0 42 ergs s~ x ) were removed. The final GOODS-N sam- 
ple consists of seven sources. Out of these, one has spec- 
troscopic redshift (z=1.86) while for the rest we adopt pho- 
tometric redshift by Le Borgne et al. 2009 (median z=1.83 
and 1.53<z<2.05). The final combined sample (LHN and 
GOODS-N) consists of 25 ULIRGs with a median z=2.01 
and with 18 out of 25 objects lying in narrow redshift range 
(1.7<z<2.3). PACS/SPIRE photometry of our sample is 
presented in Table 1 (electronic version), while IRAC, MIPS 
and mm and radio photometry is given in Table 2 (electronic 



3 DERIVATION OF FAR-IR PROPERTIES 

To derive estimates for the Lir (Ls-ioooMm) of the galaxies in 
our sample, we first convert their SED to rest-frame applying 
k-corrections and then fit the PACS and SPIRE data with 
the libraries of Chary & Elbaz (2001) (CE01) and Dale & 
Helou (2002). Results based on the two methods are in very 
close agreement indicating a median Lir=3x10 12 Lq. The 
CEOl-derived Lir for each object are summarized in Table 
1, while examples of the rest-frame SEDs along with the 
best-fit CE01 templates for two ULIRGs in our sample are 
shown in Figure 1. 

To derive the dust temperature of galaxies in our sam- 
ple, we use a single temperature modified black body fitting 
form in which the thermal dust spectrum is approximated as 
F v oc l P+f ) /( fi ( hv / kT d) _ i). This model was fit to Herschel 
data with rest-frame A > 40//m, assuming a fixed emmi- 
sivity index of /3=1.5. This wavelength cut was introduced 
to avoid fitting emission from Very Small Grains (VSGs). 
The Td of each object was obtained from the best fit model, 
based on the minimization of the \ 2 value. The uncertainty 
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for each T& value was estimated by repeating the same pro- 
cedure for random perturbations of the fitted photometric 
points within their errors (following a normal distribution). 
The best fit model for two ULIRGs in our sample are shown 
in Figure 1 (solid blue line) with the Td for each of the 
galaxies summarized in Table 1. Finally, to check whether 
the lack of PACS data for the GOODS-N sample introduces 
any systematic bias in the derived properties we repeated 
the fitting procedure for the LHN sample, excluding this 
time the PACS photometric points. The values derived with 
and without the PACS were in good agreement (< A T<j > 
= 1.9 K). 



3.1 AGN contribution to our sample 

It has been shown by previous studies that the selection 
criteria of our sample have been very successful in select- 
ing starburst over AGN dominated ULIRGs, (e.g. Farrah 
et al. 2008, Huang et al. 2009). Indeed, for five ULIRGs in 
LHN, IRS spectroscopy indicates that their mid-IR emis- 
sion is dominated by vigorous star-formation rather than 
an AGN (Fiolet et al. 2010 in prep). None of our objects 
in this field is detected by Chandra to a 0.3-2.5keV flux 
limit of 5 x 10~ 16 erg cm^s" 1 s(L^ >8.5xl0 42 erg s" 1 for 
z=2.0) (Polletta et al. 2006), while on construction of the 
GOODS-N sample all candidate objects with X-ray detec- 
tion at a flux limit of 1.95 x 10 -17 erg cm _2 s _1 were rejected 
from our analysis. As the moderate depth of the LHN X-ray 
data do not provide strong constraints on the AGN contri- 
bution we further explore this issue by the q parameter (q = 
Zo 5 [L 4 o-i2cW(3.75 x 10 12 W] -log[Li A GHz/(W Hz" 1 )], 
Helou 1985). Given that almost all of our galaxies have ra- 
dio detections, we estimate the q parameter and find a mean 
< q >=2.21 witssh intrinsic dispersion <j q =0.17. This is in 
agreement with that found by Younger et al. (2009) and 
the q of star-forming galaxies, quoted by Ivison et al (2010) 
(q — 2.40, 2cr,2=0.27). These considerations, support the con- 
clusion of Huang et al. (2009) that an AGN contributes little 
(< 10 — 20%) to the bolometric luminosity of these objects. 



4 RESULTS 

4.1 Far-IR properties and comparison with other 
ULIRG samples 

Galaxies in our sample have dust temperatures that span a 
wide range 25 ^ Td ^ 62 (K), while their luminosities vary 
by less than an order of magnitude 12.24 ^ Ioq(Lib./Lq) ^ 
12.94. The median values are T d =42.3 K, and Lir=3x10 12 
1/0, indicating a star formation rate of ~520 Mq yr" 1 (as- 
suming Salpeter IMF). It is interesting to compare these 
values to that of ULIRG samples selected by different tech- 
niques. 

We consider a large set of z-~2 SMGs (Chapman et 
al. 2005 and Kovacs et al. 2006), a sample of z~2 OFRGs 
(Casey et al. 2009) and a compilation of local/intermediate- 
z (0<z<0.98) ULIRGs (Clements et al. 2010, Farrah et al. 
2003 and Yang et al. 2007). In all these studies, the method 
to derive Td estimates is similar to ours, fitting modified 
black-body models to the far-IR photometric points. For 
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Figure 2. The Lir — Td relation for our sample (red circles). 
Included are results for local/intermediate-z ULIRGs (green filled 
triangles, Farrah et al. 2003, Clements et al. 2010, Yang et al. 
2007), high-z SMGs (blue squares, Chapman et al. 2005, Kovacs 
et al. 2006) and OFRGs (black squares, Casey et al. 2009). The 
cyan shaded area denotes the 2<r envelope of the Lir-Tj relation 
of high-z SMGs. For a given Lir, our sample span in a wide 
range of dust temperatures, bridging the "cold" high-z SMGs to 
the "warmer" local/intermediate-z ULIRGs and z~2 OFRGs. 



studies that quote Lfir instead of Lir we adopt the follow- 
ing conversion factor between the two values: Lir = 1.19x 
Lfir (Dale et al. 2001). We note that preliminary results 
by Hwang et al. 2010 (in prep) and Chapman et al. 2010 
(in prep) and Chanial et al. 2010 (in prep) indicate that the 
far-IR properties of SMGs and OFRGs when Herschel data 
are taken into account, are, in general, consistent with the 
results obtained in the pre-Herschel era. 

In Figure 2 we show the Lir — Td relation for our sample 
as compared to that of local/intermediate-z ULIRGs, SMGs 
and OFRGs. For the luminosity bin of our sample, SMGs 
have a median Td = 36 ± 8 K while OFRGs are considerably 
warmer with median Td=47 ± 3 K (Magnelli et al 2010) and 
dust temperatures similar to that of local ULIRGs. There- 
fore, it appears that the two methods select ULIRGs with 
significantly different dust temperatures, and with no signif- 
icant overlap between them. Taking advantage of the wealth 
of multi- wavelength data in GOODS-N, we find that a large 
fraction of SMGs (15/24) and OFRGs (4/5) in GOODS-N 
(Pope et al. 2006, Casey et al. 2009) that fall in the redshift 
bin (1.5<z<3.0) of our sample, satisfy the IRAC-peakers 
colour criteria, if we relax the /24 cut. 

Based on this plot there are a number of significant re- 
sults to be drawn. First of all, our observations confirm the 
existence of ULIRGs in the high-z universe with dust tem- 
perature higher than that of SMGs. Furthermore, it seems 
that the selection of high-z ULIRGs based on the detection 
of the 1.6/xm bump doen't favour a particular Td, selecting 
ULIRGs that overlap with the SMGs and OFRGs but also 
ULIRGs of intermediate Td. Indeed, we see that objects in 
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our sample range from those that are as cold as SMGs to ob- 
jects as warm as OFRGs, while a significant fraction lies in 
the intermediate region between the two samples, bridging 
the two populations. We also note that a large fraction of 
the sample falls in the Td-im relation of the local ULIRGs. 
Finally, our data indicate that the Td dispersion of high-z 
ULIRGs is larger than that of the local ULIRGs as derived 
based on IRAS observations. This discrepancy mainly arises 
due to absence of cold sources in the local universe, although 
the IRAS selection might miss existing cold sources, intro- 
ducing a bias towards warmer ULIRGs. 

4.2 SMGs: Evidence of selection bias towards 
colder ULIRGs 

There is growing evidence that ground based (sub)mm ob- 
servations introduce a systematic bias towards the detec- 
tion of cold ULIRGs. As mentioned above this was first 
discussed by Chapman et al. (2004), introducing the pop- 
ulations of OFRGs while a similar conclusion was reached 
recently by Chapin et al. (2010) using BLAST data. In Fig- 
ure 2 we showed that a fraction of IRAC peakers also tends 
to be warmer than high-z SMGs. We now ask whether these 
IRAC-peakers would be missed by the sub-mm selection. 

To investigate this, we estimate the Sgso flux densi- 
ties of our sample based on the best fit CE01 model that 
was obtained through the fitting of the Herschel photomet- 
ric points. The predicted Ssso fluxes of our sample along 
with the measured sub-mm flux of high-z SMGs are plot- 
ted over the derived Ti of the two populations in figure 
3. We also overplot tracks in constant Lir. This plot illus- 
trates that a significant fraction (~60%) of the mid-IR se- 
lected ULIRGs in our sample have Ssso flux densities lower 
than that of the SMGs, lie below the confusion limit at 
850/tm (2-3 mjy, Knudsen et al. 2008) and hence would 
be missed by ground-based (sub)mm surveys. Nevertheless, 
we also find IRAC-peakers with predicted 5*850 above the 
detection limit and which therefore should be detected in 
the sub-mm. Indeed, four of our objects in LNH (LHN1, 
LHN8, LHN16, LHN29), have been detected (S/N > 3) by 
MAMBO 1.2mm (Fiolet et al. 2009, Kovacs et al. 2010). 
For these objects we use the formula described by Ivison 
et al. (2005) to convert the observed 1.2mm to 850/mi flux 
densities and then compare these values with the predicted 
Ssso flux densities that we derived from our analysis. The 
two values are in close agreement for all objects, with a me- 
dian difference of 0.15mJy. Furthermore, we find that all 
galaxies in our sample with MAMBO observations but no 
detection (LHN0,LHN19,LHN25, Fiolet et al. 2009) have 
predicted fluxes below the detection limit. The same test 
for the GOODS-N sample reveals that our analysis, success- 
fully predicts the sub-mm fluxes of two objects with SCUBA 
850/mi detection (GN17, GN06, Borys et al. 2005, Pope et 
al. 2006). 

bf Another way to explore this issue, is a direct visu- 
alization of the SEDs of the sources. Such approach is free 
of systematics and uncertainties of SED fitting, that could 
possibly affect/bias our results. In Figure 4, we show the 
far-IR SED of four sources in LH-N that have been followed 
up by MAMBO 1.2mm observations. Two sources are not 
detected at 1.2mm while the other two have a > 3<r detec- 
tion. These sources are also chosen to have sililar Lir ( 12.5 
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Figure 3. Dust temperature versus the estimated Sg50 nu x den- 
sities of galaxies in our sample (red circles). We also include T<j 
measurements and observed S$so flux densities of high-z SMGs by 
Chapman et al, (2005) (blue squares). Solid lines represent tracks 
in constant Lir while the vertical dotted line indicates the con- 
fusion limit of current ground based submm surveys. It is evident 
that a significant fraction of our sample lies below the detection 
limit and would be missed the SCUBA-850/tm surveys, if we con- 
sider that the detection limit should be above the confusion. 



< log (Lir/Lq) < 12.7) and silimar fluxes at 250/im. It is 
evident that for the same luminosity, the far-IR SED of the 
MAMBO-undetected sources, peaks at shorter wavelengths 
compared to that of the MAMBO detected sources. This 
indicates a clear difference of the dust temperature and sub- 
sequently of the 850-1200/xm emission of the two samples, 
with MAMBO-undetected sources having warmer Td and 
considerably lower 850-1200/im flux densities. 

Recently, Kovacs et al. (2010), presented a far-IR study 
of 20 luminous z~2 mid-IR selected starbursts based on 
SHARC-2 350/tm and concluded that their properties are in- 
distinguishable from the purely SMGs population. Although 
this seems to contradict our findings this is not the case. 
Since their study focuses on IRAC-peakers with Si. 2mm > 
2mJy, their sample is biased towards the most sub-mm lu- 
minous galaxies among the mid-IR selected ULIRGs and 
hence those that are likely to share similar properties with 
the SMGs. As illustrated in Figure 3 such galaxies exist in 
our sample too. In fact, our sample shares four objects in 
common with that of Kovacs et al. (2010), for which the esti- 
mates of the far-IR properties between the two studies are in 
very good agreement. On the other hand, as we have shown 
above, due to the requirement of MAMBO detection, they 
miss a large fraction of mid-IR selected ULIRGs that have 
faint 850-1200/mi flux densities and their properties are dif- 
ferent from that of SMGs. Furthermore, although there are 
no 850-1200/im observations for some galaxies in our sam- 
ple, the fraction of galaxies with predicted Ssso above the 
detection limit is consistent with the fraction of MAMBO 
detected mid-IR ULIRGs (~40%) in the study of Lonsdale 
et al. (2009) and Fiolet et al. (2009). To summarise, Her- 
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Figure 4. The far-IR part of the SED of two MAMBO- 
detected (blue circles) and two MAMBO-undetected sources 
(black squares) from our sample. All four sources are chosen to 
have comparable Lir, and /250- The far-IR part of the SED of 
Arp220 (coral shaded area) and that of NGC2369 (grey shaded 
area) are also shown. For the MAMBO undetected soures the 
open boxes correspond to the 1.2mm flux density based on the 
SED extrapolation. The red arrow indicates the confusion limit 
of 850/im surveys. The SED of MAMBO-undetected sources 
peaks at shorter wavelenghts, indicating warmer T<j and lower 
850-1200/xm emission when compared to that of IRAM detected 
sources. This plot illustrates that among sources with compara- 
ble Lir , those with higher Tj are missed by current ground-based 
surveys. 



schel data allow us for the first time to characterize the 
far-IR properties of ~ 50% of the mid-IR selected ULIRGs 
that would be missed by ground based (sub)mm surveys 
and reveal that their properties are different from that of 
SCUBA/IRAM selected galaxies. 



5 CONCLUSIONS 

We have presented the far-IR properties of mid-IR selected 
ULIRGs at z~2 in LHN and GOODS-N fields, based on 
Herschel PACS and SPIRE observations. We showed that 
for a narrow range of luminosities, our sample spans a wide 
range of T<j, indicating that the mid-IR selection of high-z 
ULIRGs doesn't introduce a systematic bias in T&. Sources 
in our sample range from those that are as cold as high-z 
SMGs to objects as warm as OFRGs, while a significant frac- 
tion has intermediate Td, bridging the two populations. We 
also demonstrated that a significant fraction of our sample 
would be missed from (sub)mm surveys, showing that the 
sub- mm technique introduces a bias towards the detection of 
colder ULIRG sources. We confirmed the existence of star- 
forming ULIRGs at high-z that are warmer than SMGs and 
showed that the Td dispersion at high-z is larger than that 
found in the local universe. While this large dispersion in T c \ 
suggests a diversity of the physical mechanisms that drive 
the star-formation activity in the early galaxies, its origin 
remains unclear. Herschel observations of larger samples in 
the rest of the HerMES survey, will address this question 



as well as the contribution of ULIRGs to the star-formation 
density and their clustering properties. 
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Table 1. Far-IR properties of mid-IR selected z~2 ULIRGs. 



ID 


z 


Swo b 


Sieo 6 


5250 


S350 


S500 


Lir 


Td 






[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[L©1 


[K] 


LHNO 


2.01 a 


7.36 ± 3.67 


23.04 ± 4.99 


27.17 ± 1.94 


13.20 ± 1.77 


4.76 ± 2.81 


12.70 ± 0.16 


51.61 ± 2.26 


LHN1 


1.95 a 


13.37 ± 2.50 


41.78 ± 5.96 


52.51 ± 2.14 


34.21 ± 2.86 


5.02 ± 5.67 


12.99 ± 0.10 


44.43 ± 1.24 


LHN2 


2.01 


5.12 ± 1.68 


0.00 ± 0.00 


7.36 ± 1.57 


6.01 ± 1.74 


0.53 ± 1.97 


12.25 ± 0.32 


38.66 ± 4.47 


LHN3 


2.40 


8.94 ± 2.80 


0.00 ± 0.00 


15.11 ± 1.71 


15.20 ± 2.12 


14.93 ± 2.38 


12.57 ± 0.05 


33.20 ± 0.72 


LHN4 


1.72 


20.05 ± 2.83 


35.52 ± 6.37 


32.53 ± 7.87 


10.83 ± 15.04 


0.00 ± 0.00 


12.81 ± 0.11 


46.64 ± 2.30 


LHN5 


1.98 


8.70 ± 2.59 


0.00 ± 0.00 


16.86 ± 1.64 


13.12 ± 1.95 


7.07 ± 2.23 


12.52 ± 0.15 


40.00 ± 4.25 


LHN8 


2.07 a 


11.07 ± 2.53 


46.10 ± 5.23 


46.65 ± 2.03 


41.73 ± 1.67 


17.70 ± 2.19 


12.92 ± 0.10 


42.30 ± 1.04 


LHN10 


2.22 


3.01 ± 2.70 


19.29 ± 5.36 


23.60 ± 6.40 


21.09 ± 8.92 


0.00 ± 0.00 


12.86 ± 0.16 


45.74 ± 3.12 


LHN11 


2.09 


6.18 ± 2.01 


0.00 ± 0.00 


13.87 ± 2.07 


14.70 ± 2.40 


0.00 ± 0.00 


12.50 ± 0.18 


32.70 ± 0.36 


LHN16 


2.10 a 


5.47 ± 2.74 


13.44 ± 4.41 


29.35 ± 6.31 


37.53 ± 8.47 


0.00 ± 0.00 


12.57 ± 0.20 


44.11 ± 4.83 


LHN19 


1.80 


13.65 ± 2.69 


26.41 ± 5.18 


48.44 ± 2.54 


51.48 ± 3.34 


32.11 ± 4.38 


12.72 ± 0.09 


31.14 ± 0.42 


LHN20 


1.64 


6.48 ± 2.86 


13.90 ± 4.39 


13.16 ± 1.63 


5.65 ± 1.68 


0.00 ± 0.00 


12.35 ± 0.24 


50.81 ± 4.78 


LHN24 


2.15 


2.82 ± 4.05 


21.95 ± 6.87 


16.57 ± 1.77 


15.74 ± 2.11 


12.31 ± 3.81 


12.68 ± 0.15 


38.64 ± 3.64 


LHN25 


2.18 


4.55 ± 2.42 


13.87 ± 4.21 


14.60 ± 3.38 


3.42 ± 3.83 


9.82 ± 2.13 


12.56 ± 0.26 


38.82 ± 3.61 


LHN27 


2.23 


5.23 ± 3.44 


18.90 ± 4.68 


15.47 ± 2.44 


10.34 ± 2.39 


12.49 s ± 3.57 


12.61 ± 0.28 


52.00 ± 0.68 


LHN29 


1.96 a 


0.00 ± 0.00 


12.83 ± 4.11 


27.47 ± 1.75 


32.45 ± 2.74 


18.45 ± 3.80 


12.51 ± 0.19 


30.51 ± 1.34 


LHN30 


1.56 


9.53 ± 2.82 


22.48 ± 5.10 


20.80 ± 2.00 


11.37 ± 1.62 


0.00 ± 0.00 


12.43 ± 0.18 


45.75 ± 2.43 


LHN31 


3.03 


5.38 ± 2.79 


12.77 ± 4.22 


21.98 ± 1.61 


14.76 ± 2.51 


0.00 ± 0.0 


12.94 ± 0.17 


62.10 ± 5.87 


GN18 


2.05 


- 


- 


9.46 ± 1.09 


11.72 ± 1.62 


0.00 ± 0.00 


12.38 ± 0.14 


30.46 ± 2.44 


GN32 


1.83 


- 


- 


41.36 ± 1.02 


41.46 ± 1.62 


0.00 ± 0.00 


12.72 ± 0.14 


35.70 ± 2.86 


GN34 


1.83 


- 


- 


15.12 ± 1.37 


13.64 ± 2.10 


0.00 ± 0.00 


12.28 ± 0.14 


35.44 ± 2.83 


GN35 


1.83 


- 


- 


37.47 ± 2.04 


32.22 ± 4.45 


0.00 ± 0.00 


12.69 ± 0.14 


39.82 ± 3.19 


GN44 


1.74 


- 


- 


13.51 ± 1.55 


16.18 ± 2.65 


0.00 ± 0.00 


12.24 ± 0.13 


34.46 ± 2.76 


GN46 


1.66 


- 


- 


41.45 ± 1.45 


52.20 ± 2.06 


26.40 ± 4.24 


12.71 ± 0.14 


36.13 ± 2.89 


GN58 


1.52 


- 


- 


23.59 ± 1.07 


12.59 ± 2.05 


0.00 ± 0.00 


12.41 ± 0.14 


44.58 ± 3.57 



a IRS spectroscopy by Fiolet et al. (2010) 

b No available PACS data for GOODS-N in this study 



Table 2. Summary of ancillary data. 
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ID 


RA 


DEC 


<?3.6 


S4.5 


S5. 8 


S&.o 


S24 


SlAGHz 








I/^Jy] 


[^Jy] 


t/^Jy] 


[j«Jy] 


0*Jy] 


t/^Jy] 


LHNO 


161.127548 


58.921799 


41.3 ± 0.8 


49.5 ± 1.1 


55.3 ± 4.0 


53.9 ± 3.8 


781 ± 24.0 


101.7 ± 14.2 


LHN1 


161.487091 


58.888611 


26.4 ± 0.7 


33.6 ± 1.0 


39.5 ± 3.6 


36.9 ± 3.6 


684 ± 24.1 


314.8 ± 19.1 


LHN2 


161.376022 


58.920658 


28.4 ± 0.7 


34.2 ± 0.7 


40.2 ± 3.5 


27.1 ± 2.6 


375 ± 22.3 


29.0 ± 7.5 


LHN3 


161.415726 


58.906940 


30.9 ± 0.6 


39.5 ± 0.9 


47.2 ± 3.2 


47.8 ± 3.0 


485 ± 24.1 


46.9 ± 4.3 


LHN4 


161.545685 


58.879189 


52.9 ± 0.9 


66.9 ± 1.1 


57.0 ± 3.6 


52.1 ± 3.4 


401 ± 25.0 


160.2 ± 10.7 


LHN5 


161.160263 


59.075150 


32.3 ± 0.6 


37.8 ± 0.9 


30.0 ± 3.1 


33.5 ± 3.4 


375 ± 23.3 


51.3 ± 9.8 


LHN8 


161.661163 


58.936852 


29.5 ± 0.7 


38.8 ± 1.0 


58.2 ± 3.7 


42.0 ± 3.5 


662 ± 23.4 


159.5 ± 9.9 


LHN10 


161.525223 


59.141270 


45.4 ± 0.7 


58.0 ± 1.0 


60.2 ± 3.0 


45.7 ± 3.2 


194 ± 25.6 


238.9 ± 16.2 


LHN11 


161.231583 


59.167912 


48.9 ± 0.8 


59.0 ± 1.1 


54.4 ± 3.3 


36.0 ± 3.3 


258 ± 24.8 


45.1 ± 10.0 


LHN16 


161.824844 


59.042171 


35.8 ± 0.8 


50.0 ± 1.1 


77.4 ± 4.0 


52.1 ± 3.8 


567 ± 26.4 


55.1 ± 5.4 


LHN19 


161.507462 


59.154690 


72.1 ± 1.0 


91.3 ± 1.5 


93.5 ± 3.9 


74.0 ± 3.9 


1011 ± 23.1 


82.2 ± 5.7 


LHN20 


161.676163 


59.069839 


53.2 ± 0.8 


61.6 ± 1.2 


57.3 ± 3.3 


42.9 ± 3.5 


307 ± 23.9 


70.6 ± 5.0 


LHN24 


161.843964 


59.019981 


31.9 ± 0.6 


41.0 ± 1.0 


41.1 ± 3.1 


38.0 ± 3.3 


439 ± 23.3 


73.7 ± 9.4 


LHN25 


161.933746 


59.106892 


36.3 ± 0.5 


43.5 ± 0.8 


46.8 ± 2.8 


36.9 ± 3.3 


656 ± 25.7 


54.9 ± 13.7 


LHN27 


161.729813 


59.191101 


28.9 ± 0.6 


36.8 ± 0.7 


40.0 ± 3.5 


38.4 ± 3.0 


337 ± 22.5 


36.4 ± 10.5 


LHN29 


161.909683 


59.169449 


44.4 ± 0.6 


51.0 ± 0.8 


52.3 ± 3.1 


47.7 ± 3.3 


688 ± 24.0 


69.2 ± 9.3 


LHN30 


161.944641 


59.257740 


50.5 ± 0.6 


59.4 ± 0.9 


40.0 ± 2.9 


50.3 ± 3.3 


341 ± 23.4 


87.6 ± 16.8 


LHN31 


161.935730 


59.210369 


35.8 ± 0.6 


44.5 ± 0.9 


39.5 ± 3.1 


35.8 ± 3.4 


404 ± 21.8 


0.0 ± 0.0 


GN18 


189.262684 


62.142494 


11.4 ± 0.1 


13.9 ± 0.1 


10.3 ± 0.5 


13.3 ± 0.6 


230 ± 7.2 


- 


GN32 


189.256739 


62.196195 


53.9 ± 0.1 


70.1 ± 0.1 


8.8 ± 0.3 


57.9 ± 0.4 


716 ± 8.1 


- 


GN34 


189.399541 


62.345261 


29.3 ± 0.1 


37.5 ± 0.1 


25.3 ± 0.5 


28.2 ± 0.5 


178 ± 5.4 


- 


GN35 


189.076657 


62.264067 


14.7 ± 0.1 


19.6 ± 0.1 


22.2 ± 0.4 


20.4 ± 0.4 


314 ± 5.3 


- 


GN44 


189.074144 


62.235591 


49.2 ± 0.1 


55.6 ± 0.1 


35.9 ± 0.4 


40.2 ± 0.5 


428 ± 7.2 


- 


GN46 


189.297273 


62.225206 


37.9 ± 0.1 


45.0 ± 0.1 


14.8 ± 0.3 


37.8 ± 0.4 


534 ± 8.6 


- 


GN58 


189.294242 


62.376245 


38.7 ± 0.1 


47.9 ± 0.1 


12.4 ± 0.7 


46.4 ± 0.6 


383 ± 4.6 


- 
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Table 3. SHARC-350^m and MAMBO 1.2mm flux dcnities of our LHN sample 



ID 


S(350fim) a 


S(1.2mm) b 




[mjy] 


[mjy] 


LHN1 


39.7±5.9 


3.08±0.58 


LHN8 


31.9±4.9 


2.13±0.71 


LHN16 


49.7±6.5 


2.66±0.78 


LHN29 


31.8±5.9 


2.48±0.74 



a Flux densities by Kovacs et al. (2010) 
b Flux demsities by Fiolet et al. 2009 



